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Abstract -We prove existence of discrete solitons in infinite parity-time (PT-) symmetric lattices 
by means of analytical continuation from the anticontinuum limit. The energy balance between 
dissipation and gain implies that in the anticontinuum limit the solitons are constructed from 
elementary PT-symmetric blocks such as dimers, quadrimers, or more general oligomers. We 
consider in detail a chain of coupled dimers, analyze bifurcations of discrete solitons from the 
anticontinuum limit and show that the solitons are stable in a sufficiently large region of the 
lattice parameters. The generalization of the approach is illustrated on two examples of networks 
of quadrimers, for which stable discrete solitons are also found. 
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Introduction. — Energy localization in lattices is a 
fundamental topic. It received its particular significance 
after the prediction of the intrinsic localized modes [1] 
and the subsequent rigorous proof of the existence of such 
modes [2] by analytical continuation from the anticontin- 
uum limit when the coupling of the nearest neighbors is 
weak. Nowadays the topic is very well elaborated and 
numerous physical applications, including nonlinear op- 
tics and [3] Bose-Einstein condensates [4] have been thor- 
oughly studied. One of the most popular model appear- 
ing in description of these physical phenomena, which is 
also a widely accepted testbed for mathematical analy- 
sis of the anticoninuum limit, is the discrete nonlinear 
Schrodingcr equation (DNLS) [5], also known as the dis- 
crete self-trapping equation [6]. 

More recently, particular attention was paid to the 
DNLS with gain and losses. Such models naturally appear 
in the optical context of arrays of amplifying and absorb- 
ing waveguides [7]. If gain and losses are adjusted to create 
the refractive index profile having symmetric real and anti- 
symmetric imaginary parts [8], such systems have parity- 
time (PT) symmetry and may have pure real spectrum. 
Originally the idea about existence of pure real spectrum 
of complex potentials obeying the PT-symmctry was in- 
troduced in [9] questioning the fundamentals of the quan- 
tum mechanics. It turned out however that the most direct 
applications of the PT symmetry today can be found in 
the discrete optics. Namely, in such systems, and more 



specifically in two coupled waveguides (one with dissipa- 
tion and another with gain) the phenomenon was imple- 
mented experimentally |10j . 

Many detailed studies of PT-symmetric DNLS equa- 
tion were already developed for lattices with a finite 
number of sites. In particular, there were considered 
periodic oscillations in a system of two oscillators (a 
dimer) [llj ; stationary nonlinear modes for four oscilla- 
tors (a quadrimer) |12|ll3j: the relation between one- and 
two-dimensional finite "PT-symmetric networks [13j . and 
detailed analysis of two-dimensional plaquettes [TJ]. The 
transition to the limit of an infinite number of sites was 
investigated in [15]. It was shown that in this limit the 
PT-symmetry breaking occurs at gain-loss coefficient ap- 
proaching zero. Stable discrete solitons in the infinite 
chain of PT-symmetric DNLS oscillators with alternating 
coupling were discovered numerically in |16j . However, the 
solitons obtained in |16j , displayed oscillations and neither 
analytical proof of the existence nor the number of possi- 
ble families of solutions were clarified, so far. 

In this Letter, we give an analytical proof of the exis- 
tence of discrete solitons in PT-symmetric DNLS lattices 
with alternating coupling coefficients and classify differ- 
ent solution families and their stability near the anticon- 
tinuum limit. In particular, we report multistability of 
localized modes, that is, the existence of two or more sta- 
ble solutions with the same energy and the same lattice 
parameters (but having different shapes). 
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Fig. 1: Schematic presentation of the PT-symmetric array of 
waveguides with gain, "+", and losses, "— ", and with the no- 
tations used in the text. 

The DNLS equation with alternating coefficients of gain 
and loss can be viewed is a discrete (tight-binding) limit 
of a continuous 'PT-symmetric lattice. Stable solitons in 
such systems have been found in the presence of only lin- 
ear [17], only nonlinear [18], and both linear and non- 
linear [19] PT-lattices. The solutions considered in this 
Letter can be viewed as discrete counterparts of the men- 
tioned solitons. 

We consider the "PT-symmetric DNLS equation 

i <kn + _ qn j _|_ Cn+1 (g n _ 1 - q n ) - g\q n \ 2 q n 

at 

+i{-l) n+1 iqn = 0, (1) 

where the positive constants c„ = kq for n = 2p and c n = 
Ki for n = 2p + 1 describe the two alternating couplings 
(kq and «i, with fto,i > 0) between neighbor sites, and 
it is assumed that all odd (even) sites have loss (gain) 
which is characterized by factor 7 > [see fig. Q]. In 
the context of optical applications our model describes an 
array of waveguides with gain and losses. Then q n is a 
dimensionless field in the waveguide n, and t means the 
propagation coordinate. 

Let us first briefly address the most important features 
of the underlying linear problem, which can be formally 
obtained from eq. §Q by setting g = 0. Identifying solu- 
tions of the linear problem in the form of Floquet-Bloch 
modes (q 2n , 9271+1) = (u,v)e tkn ^ Ko+K - 1+ ^ t , one can re- 
cover that PT symmetry is unbroken on the infinite lattice 
if [IB] 

|«o - «i| > 7- ( 2 ) 

In this case for any real k the corresponding eigenvalue /x 
is real. More precisely /j, 2 lies in the interval 

(k -ki) 2 -7 2 < M 2 < (k + ki) 2 -7 2 , (3) 

i.e. fi belongs either to a positive or to a negative spec- 
trum band. If inequality ([2]) does not hold, we say that 
PT symmetry is broken as there exist eigenvalues [i with 
nonzero imaginary parts. 

Returning to the nonlinear problem, without loss of gen- 
erality we impose g = ±1. Furthermore, by analogy with 
the conservative DNLS equation (see e.g. [20]), one can 
verify that there exists the symmetry reduction as follows. 
If q n is a solution of © for g = 1, then (-l) n q n e- 2 ^ K «+ K ^ t 
is a solution of ([TJ for g = —1. This reduction allows us 
to restrict further consideration to the case g — 1 only. 

Anticontinuum limit. — We are concerned with the 
stationary solutions, i.e. solutions whose dependence on 



time is given by q n (t) ~ e~ lEt , where E is a constant which 
is termed an energy (or a propagation constant in optical 
applications). Then, in the case of the conservative DNLS 
equation, the anticontinuum limit corresponds cither to a 
coupling between the two neighbor sites tending to zero 
or to the energy tending to the infinity, the two limits 
being equivalent, i.e. mapped to each other by simple 
transformation (see e.g. [20]). The same is true for the 
model ([1]), however with two important changes. 

First, in the presence of the dissipative term i(— l) n "iq n 
the elementary cell of the DNLS equation (JIJ is composed 
of two sites, one with gain and the other one with loss 
(even if kq = k\). Therefore, the anticontinuum limit 
must be formulated in terms of dimers, rather than single 
sites. Thus, the anticontinuum limits can be identified 
as either kq or «i to be small enough. Without loss of 
generality, we fix kq = 1. Then the anticontinuum limit 
corresponds to Hi = 0. 

In physical applications, both coupling constants kq and 
K± are usually fixed. Then the anticontinuum limit can 
be realized at the limit of large energy E. However, and 
this is the second distinction from the conservative DNLS 
equation, although the small parameter can be obtained 
in the limit of large energy, E cannot be scaled out from 
the main equations if 7 ^= 0. 

Using the above considerations, it is convenient to 
rewrite the main model ([1]) in terms of variables 

( Q2n(t) \ _ ( u n \ 4 ( Ko+Kl+AI ) t (A] 
\<hn+l(t) )~\Vn) e ' W 

where [i is a constant and we assume that u n and v n do 
not depend on t and satisfy zero boundary conditions: 
u n ,v n — > as ri — > ±00. Then for kq = 1, k\ = e, 
and g = 1, the main model can be rewritten in the ma- 
trix form for w n = (u n ,v n ) T ("T" stands for the matrix 
transposition) : 

Hw n + e (o-_u>„_i + (T + w n+1 ) = F(w n )w n , (5) 

where 

H = f M - *7 1 \ 
V 1 A* + *7 J ' 

ff - = (o J)' fT+ = (! 0)' 

and F(w n ) = diag(|u„| 2 , K| 2 ). 

Single-dimer state. — First, we address the simplest 
case when in the anticontinuum limit e = 0, only one 
central dimer is excited, i.e. Wq ^ 0, whereas w n = 
for n/0. Then for any n ^ eq. ([5]) is automatically 
satisfied. For the central dimer, i.e. at n = 0, we assume 
that the following PT symmetry reduction holds: uq = vq 
and arrive at the following equation |11| : 

(H - ij)u + u = K| 2 uo- (6) 
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The latter equation has two branches of solutions, which 
exist for all 7 G [0, 1): 

ut=A ± e i *±, (7) 

where A\ = fx±^l — 7 2 , sin(2</?±) = —7, and cos(2<^±) = 
— 7 2 . The branches exists for p > /i± where 
/i± = ±^1 — 7 2 . The physical difference between 
these branches becomes evident if we introduce the lin- 
ear momentum p n = 2Rc(q n q n +i) and the current j n = 
2Im(g„<7„ + i) per unit cell; respectively p = J2 n P n are 
j = Y] j n are the total momentum and current carried 
out by the solution. Then, branch (branch Uq) corre- 
sponds to the linear momentum and current, between the 
two sites of the dimer, having the same (opposite) direc- 
tions. Note that eq. (O coincides with the equation for 
the stationary solutions of the parametrically driven NLS 
equation [2T[|2"2]. 

Looking for continuation of the solution ([7]) from the 
anticontinuum limit (i.e. from e = to e > 0) it is nat- 
ural to suppose that for e > the solitons also obey the 
symmetry reduction on the whole infinite lattice, i.e. 

U n = V- n , V n =U- n , H = 0, ±1, ±2, (8) 

It allows one to restrict the consideration only to n > 
0. At the central dimer, Wq, one can introduce the real 
coordinates (ao,6o) such that uq = vq = ao + ibo and 
rewrite §5§ for n = as follows 

f na + 760 + ao + eRe(iti) = (a 2 , + 6g)a , 
\ M 6 o - 7 a o - b - elm(ui) = (a 2 , + bl)b . 

For n > 1, we still use the complex- valued coordinates: 

f (/i - ij)u n + V n + ev n -l = \u n \ 2 U n , 
\ (n + ij)Vn+Un+eUn + l = \Vn\ 2 V n . 

Now the system (|5j)- (flTj|) is smooth with respect to pa- 
rameter e and the solution vector. At e = 0, we have the 
limiting solution: u n — v n — for all n = 1,2,..., while 
ao and &o are given by one of the two possible solutions (J7|) 
for 7 £ [0, 1) and fi > fi±. To apply the implicit function 
theorem arguments, we need to show that the Jacobian 
operator of the system ([<l))- (fTU| at e = is invertible at 
the limiting solution. Furthermore, the solution can be 
analytically continued from the anticontinuum limit until 
a critical value e cr > for which the Jacobian operator 
becomes non-invertible. 

In the case of the conservative DNLS equation (7 = 0) 
rigorous estimates for e cr can be obtained analytically [21 
120) . Due to mathematical constrains such estimates are 
typically lower than the practically achievable values of e cr 
for which the localized solutions exist, on the one hand, 
and on the other hand require more elaborated analytical 
study. Therefore, here we restrict the consideration only 
to the proof that the analytical continuation is possible 
and study the continuation numerically. 



For e = 0, the lattice consists of a set of uncoupled 
dimers. For n = 1,2, . . ., the limiting Jacobian operator 
of the system (fTO)) is nothing but H and thus is invertible 
for /i ^ fi± since det(H) = /i 2 + 7 2 — 1. 

At the central dimer n = 0, the limiting Jacobian oper- 
ator of the system © is given by the 2x2 matrix: 

t _( -2al-jb Q /a 7 - 2a b Q \ , 
° V -7-2a 6 -2bt + 1 a /b J ' {LL) 

where eq. (j6]) has been used. The matrix Jo is invert- 
ible if ao&o 7^ and a 2 , ^ 6q. This gives the constraints 
A± 7^ and cos(2ip±) ^ in the limiting solution (JT)), or 
equivalently, A\ ^ {0,/x}. The constraints arc satisfied 
for any 7 6 [0, 1) and /1 7^ [i±. Hence, for any 7 and \x 
that satisfy the found invertibility conditions, solutions 
of the dimer problem give birth to two branches of local- 
ized discrete solitons on the infinite T-'T'-symmetric lattice. 
These branches, which will be respectively denoted as 
[sec fig. El below], are parameterized by e, they persist at 
least for all small e, and for small e the solitons from the 
branches T 1 -^ are nearly localized at the central dimer u;o- 

Multi-dimer states. — One can also consider the 
case when the solution in the anticontinuum limit con- 
sists of several excited dimers. Say, for the case of two 
dimers, one can consider branches r( + ' + ) or r( - ' - ), which 
at e = correspond to two dimers occupying two con- 
secutive central cells n = and n = 1. More complex 
configurations consisting of N excited dimers can also be 
continued from the anticontinuum limit. Even more gen- 
erally, there exist branches like r^ + ' 0,+ ', r^ - ' ' - ), which in 
the anticontinuum limit correspond to two dimers placed 
at n = — 1 and n = 1 separated with an "empty" dimer at 
n = (i.e. Wq = for e = 0). However, the continuation 
is not possible for arbitrary choice of N central dimers. If 
we consider the sequence a = (ai, a% . . . , «jv) consisting 
of N symbols ai,...^ £ {+,—,0}, then existence of the 
branch T a is only possible provided that a p = aN+i-p 
for p = 1,...,N. The latter requirement ensures that the 
configuration, which is chosen to be continued from e = 
to e > 0, obeys a VT symmetry reduction ©. For ex- 
ample, for N = 3, the family r( - ' +,+ ' can not be con- 
tinued from the anticontinuum limit. However the family 
r(~' + '~) does persist for small e. Invertibility of the Ja- 
cobian operators corresponding to such multi-dimer solu- 
tions can be proven using similar technique as for the case 
N = 1 considered above. We note that the relevance of 
the symmetry properties for possibilities of continuation 
of the branches of soliton solutions is similar to the case 
of parametricalally driven NLS systems [22] . 

Stability in the anticontinuum limit. We shall 
also address linear stability of solitons belonging to the 
branches bifurcating from the one-dimer states in the 
anticontinuum limit. For e = and n 7^ 0, the dimers are 
decoupled and the stability of the zero solution is deter- 
mined by the spectrum of the matrix H which has real 



p-3 



V. V. Konotop et al. 



eigenvalues if 7 < 1. Hence the zero solution for n ^ is 
stable. Passing from e = to e > 0, eigenvalues A group 
together into bands of continuous spectrum. For small 
positive e, these bands are situated in the neutrally stable 
imaginary axis, and they are separated from each other 
and from zero if 7 G [0, 1), fx 7^ f/,±, and fi ^ 0. 

Thus, to ensure stability of the single-dimer soliton, we 
have to address only the stability of the central dimer wo. 
Considering a perturbed solution w + i(>o eXt + i/'i e ' U and 
linearizing the equation with respect to i/>o,i, for e = we 
obtain the eigenvalue problem 



Lq Li 



4>i 



(12) 



where L = 2 diag(|ug | 2 , \uf | 2 ) - H and L\ = 
diag((up ) 2 , (Sq ) 2 ), where are defined by eq. for 
branches rW, Because the nonlinear system ([5]) admits 
gauge invariance, A = is a double eigenvalue of the eigen- 
value problem ([12")) . As a result, the characteristic poly- 
nomial -D(A) can be factorized by A 2 and reads as (see 

also |H3): D(X) = A 2 (a 2 + 8^T~^(^T~^ ± M /2)) , 

where eq. ([7]) has been used. This expression shows that 
for e = and f <1 the solitons from branch r( + ) are sta- 
ble for any (U > /x_ and 7 G [0, 1). Solitons of the branch 
r' - ) are stable for e = and e <C 1 only if ll+ < fi < 2[i + 
and unstable with a positive eigenvalue for fi > 2/z+. 

Numerical results. — Turning now to the numerical 
study of the discrete solitons in the infinite lattice, we have 
computed bifurcations of families from the anticon- 
tinuum limit e = 0, considered their continuations to the 
domain e > 0, and examined stability of the found solitons. 
The results are conveniently visualized in the plane (P, e) , 
where P = J2 n (\ Un \ 2 + l w «| 2 )> which in optics corresponds 
to the total energy flow. In fig. [5] we show the results for 
different fi and 7. We recall that the branch r' + ) (r'~') 
is found by means of continuation starting from the dimer 
solution Uq (uq ) given by eq. ([7]). We tested several values 
of fi and 7 and in all cases numerical results for stability 
of the families for small e were in agreement with the 
above linear stability analysis. For example, branch 
is stable for fi = 1 and 7 = 0.1, but is unstable for any 
other considered values of fi and 7 on fig. [2j 

At a certain value of e = eo, the norm of the solitons 
belonging to the branch r( - ) vanishes, i.e. P — > 0. Since 
this is the linear limit, at the point P = the parameters 
obey the relation /i 2 = (1 + eo) 2 — 7 2 , which means that 
the solution branch ends up at (or bifurcate from) the 
edge of the linear spectrum [see eq. ([3])]. Respectively, 

eo = vV +7 2 - I- 

Bifurcation of the discrete solitons from the edge of the 
linear spectrum becomes particularly evident if we employ 
representation on the plane (P, fx), which is to be obtained 
for fixed 7 and e. Then, as shown in fig. [31 the found 
discrete solitons constitute continuous families, which is a 
frequent feature of nonlinear 'PT-symmetric systems 113] . 
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Fig. 2: P vs. e. Panels in the top (bottom) rows correspond to 
/i = 1 (^t = 10). Panels in the left (right) column correspond to 
7 = 0.1 (7 = 0.9). Stable (unstable) solitons are shown by solid 
blue (dotted red) lines. Notice that panel (a) has logarithmic 
scale in the horizontal axis and broken vertical axis. Vertical 
dotted line in panel (c) corresponds to e= 0.8; see also fig. [3] 



We notice that in the context of a parametrically driven 
NLS system, the connection of the soliton branch with the 
continuous spectrum was reported in |21j . 
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Fig. 3: P vs fi for 7 = 0.1 and different e. Vertical shadowed 
domains show the bands of the linear spectrum. Vertical dotted 
line in the left panel corresponds to fj, — 10; see also fig. [2] 

The branch r^ + ' is stable for sufficiently small e for all 
considered [i and 7 (in agreement with the linear stabil- 
ity analysis). For fj, = 10 [see fig. |2(c)-(d)] solitons of the 
family r^ + ' lose stability (the part of the continuous spec- 
trum leaves the imaginary axis and becomes unstable) at 
e = 1 — 7, i.e. at the VP symmetry breaking bifurcation 
[see eq. At e = 1 + 7, the VT symmetry is restored 
[see eq. fl2J], but the branch r' + ' does not become stable, 
because isolated unstable eigenvalues persist in the spec- 
trum of linearization. Altogether, the branch r^ + ) displays 
"snaking" behavior and has several alternating domains of 
stability and instability. Finally, branch T^' returns to 
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Fig. 4: Amplitude and currents for the unstable soliton at e = 
1, P « 22 [panel (a)] and for stable soliton at e — 3.32, P ~ 56 
[panel (b)]. For both shown solitons fj, = 10 and 7 = 0.1. 
Filled and empty circles correspond to sites with gain (i.e. u„) 
and losses (i.e. v n ), respectively. Arrows show directions and 
amplitudes of the largest currents in the system. 
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Fig. 5: Two examples of PT-symmetric networks, which con- 
sists of a set of disconnected quadrimers in the anticontinuum 
limit. 



the anticontinuum limit by means of coalescing with the 
branch r©' +,_ ) bifurcating from the three-dimer state in 
the anticontinuum limit. 

For large fi, e.g. \i — 10 in fig. [5Jc)-(d), branches 
can be continued into the region e£ (1 — 7, 1 + 7), where 
PT symmetry is broken. In particular, solitons exist at 
e = 1, i.e. kq = Ki, [fig. Ufa)]; such solitons, however, are 
unstable. A stable soliton is shown in fig. EJb). 

Generalizations. — The developed approach can be 
applied to the case when the elementary cell of a network 
is a more complex PT-symmetric cluster (than the dimer) . 
To illustrate this, we now briefly address the anticontin- 
uum limit for two networks of quadrimers, i.e. clusters of 
four sites w n = (w„ ,Wn ,Wn ,Wn) T , whose examples 
are shown in fig. [5] To describe the network on fig. EJa), 
we can still adopt cq. ([5]), where 



H 



( n-vy 
1 


V 



1 

fj, - i-f 
1 




\ 





1 

M + *7 1 

1 fi + i"f j 



(13) 



the nonlinearity is given by 



F(w n ) = diag 



,(!)| 2 l,„( 2 )| 2 U„(3)|2 |,„(4)|2 



and <r± are now 4x4 matrices whose only nonzero elements 
are (cr_)i4 = (er+)4i = e. The matrix H is invcrtiblc unless 

M 2 = §-7 2 ±iV5- 167 2 - 

In the anticontinuum limit, defined by e = 0, the net- 
work shown in fig. [5] (a) consist of a set of disconnected 
"PT-symmetric quadrimers. Here we consider the simplest 
case, when at e = only one central quadrimer is excited, 
i.e. w n — for n 7^ 0, and look for continuation of this 
solution to e > 0. 

To prove the possibility of analytical continuation as 
above, we concentrate on 'PT-symmetric solutions (i.e. 



obeying the symmetry Wn = and wffl — 

This allows us to restrict the consideration to the semi- 
infinite matrices with n > 0. Moreover the invertibility of 
H for /i 2 7^ I — 7 2 ± i y5 — I67 2 ensures the continuation 



provided the Jacobian matrix for the central quadrimer is 
invertiblc. 

For n = and e = 0, the central quadrimer obeys the 
system of four equations [see eq. ([5])] 



Hw = F(w )w , 



under the symmetry: Wq 1 ' = and wjf 
the nonlinear system (fl"4| generally does not admit an 
explicit analytical solution [in contrast to the dimer case 
([7])]. its properties are well studied. In particular, families 
of its nonlinear modes, bifurcation diagrams, and some 
exact solutions have been reported [T2lll4j . 

The network in fig.[5ja) is characterized by two types of 
the PT symmetry, the local and global ones. We say that 
the lattice is locally PT-symmetric if the system (fH)) is 
PT-symmetric in the limit e = 0. On the other hand, 
we say that the lattice is globally PT-symmetric if the 
infinite network ([5]) with the matrix H in f| 13[) is PT- 
symmetric for e ^ 0. The network in fig. [SJa) consists 
of quadrimers which have unbroken local PT symmetry, 
at least for small 7. For e > the infinite system has 
unbroken global PT symmetry allowing for stable discrete 
solitons. An example of a stable discrete soliton for this 
network is shown on fig.^a). 

We shall now consider the network, which consists of 
clusters whose local PT symmetry is broken. However, 
proper choice of the coupling e > makes the infinite 
network to possess unbroken global PT symmetry. An 
example of such network is presented in fig. [5jb). For this 
network, we can still work with cq. ([5]) , where 



(14) 



r,(3) 



u>o ' ■ While 



H = 



( fj, — i-y 


1/2 
V 1/2 





-1/2 
-1/2 



\ 



1/2 1/2 

-1/2 -1/2 
/i + «7 

fi + 17 J 



(15) 



Local PT symmetry is broken for any 7 because eigenval- 
ues of H are complex for any 7 > 0. (We notice that the 
local PT symmetry can be fixed if we add the diagonal 
matrix diag(l, —1, 1, —1) to H . In this case, both networks 
shown in fig. [5] have equal spectra |13j). 
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Operator H is invertible unless /i = ± 1/1/2 — 7 2 or 
/i = 7 = 0. Existence of analytical continuation of the one- 
quadrimer state from e = can be shown using the same 
ideas as the presented above. The only essential difference 
is that now PT-symmetric reduction is as follows: at the 
central quadrimer, we set w { 



(i) 



-w 



(4) 



while for n ^ 0, we set Wn = -wl n ana Wn' = W! „ 

Because local VT symmetry is broken for any 7, the 
global VT symmetry of the infinite network is also broken 
for small e. Therefore, all the solitons bifurcating from the 
anticontinuum limit are unstable at least for sufficiently 
small e. However, by increasing the coupling parameter 
e, the global VT symmetry is restored and the network 
in fig. GJb) may possess stable solitons. An example of a 
stable discrete soliton for this network is shown in fig.EJb). 
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Fig. 6: (a) Stable soliton for the network in fig. EJa) at e = 0.5, 
7 = 0.25 and /j, = 2. (b) Stable soliton for the network in 
fig- Ob) at e = 1.6, 7 = 0.1 and [i = 10. Filled (empty) circles 
correspond to sites with gain (losses). 



Conclusion. — In this Letter we have shown that the 
idea of analytical continuation from the anticontinuum 
limit can be extended to the networks of PT-symmetric 
clusters, offering analytical proof of the existence of lo- 
calized discrete solitons. Such solitons obey the PT- 
symmetric shape and can be found stable. As particu- 
lar examples, we considered in details the chains of PT- 
symmetric dimers and the networks of PT-symmetric 
quadrimcrs. 

The considered systems allow for further straightfor- 
ward generalizations, say to chains of clusters where there 
exist more than one link among the neighbor ones, like the 
chain of dimers with pairwise coupling considered in [23] 



or the chain of oligomers, i.e. clusters with more than four 
sites. Furthermore, the approach of continuation from the 
anticontinuum limit can be used for developing of a clas- 
sification of intrinsic localized modes, as well as analytical 
theory of the nonlinear stability of such modes. 
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